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trie wave in a stratified medium and the quantum mechanical Schrodinger equation in 
a piece-wise potential, we investigate the analog of the so called particle limit in quan- 
O" 1 turn mechanics. It is shown that in this limit the resonance phenomena are lost since 

individual reflection and transmission terms no longer overlap. The result is a stationary 
zebra-like response with the intensity in each stripe calculable. 
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(N : I. INTRODUCTION 

, There exist many unanswered questions in potential theory quantum mechanics. Amongst these is 

' the existence of multiple diffusion phenomena [1,2], the Hartman effect with its apparent violation 

\ of causality [3,4], the importance, if any, of wave packets in oscillation phenomena [5,6]. Most of 

these lack direct experimental measurements. It is therefore extremely instructive to study that class 
of Maxwell equations which are analogous to the Schrodinger equation. These analogies are well 
known in optics e.g. we often find references to tunneling and resonance phenomena. However, their 
5^ ' relevance to quantum mechanics has not been fully exploited. This paper studies an example of the 

above analogy. 

From the Maxwell equations, we can obtain differential equations which the electric and the 
magnetic vector must separately satisfy [7]. For example, for the electric field E, in the case of no 
charges or currents, one has 

v'-E- ^ d tt E+ (Vlnu) x (V x E) + W(E ■ Vine) = . (1) 

c 2 

The corresponding equation for the magnetic field H is obtained by making the changes e fi 
and E — * H . We shall confine our attention to the study of a medium characterized by a real (no 
attenuation) refractive index whose properties are constant throughout each plane perpendicular to 
the chosen direction z, stratified medium [8], 

n(z) = { n, for z < , n u for < z < L , n UI = n 1 for z > L } , (2) 

and for which [i assumes the same value in all three regions. By taking the plane of incidence to be 
the y-z plane, for a monochromatic, exp(— iuot), transverse electric wave (E 2 , 3 = 0), Eq.(fT]) reduces to 

d yy E 1 {y,z) + d zz E 1 (y,z)+n 2 {z)k 2 E 1 {y 1 z) = , (3) 

with k = cu/c and where we have used the fact that V • {eE) = implies that E 1 is a function of y and 
z only. The components of the magnetic vector can be determined by using V x E — — dt (fiH)/c, 

{H 2 (y,z),H 3 (y,z)}^^-{-d z E 1 (y,z),d y E 1 (y,z)} . (4) 
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The geometry of our problem is schematically represented in the following picture: 




z 



L 



Looking for a separable solution of the form 

E 1 (y, z) = U(z) exp(? n, sin 9 ky) , 

with 9 representing the incidence angle, we obtain the following second-order linear differential equa- 
tion for U(z), 

U"(z) + [n 2 (z)-n 2 +n 2 cos 2 6} k 2 U(z) = . (5) 

This equation is formally identical to the one-dimensional Schrodinger equation when the factor 
which multiplies U(z) is replaced by 2m(E — V)/h . As for Schrodinger the solutions are oscillatory 
(travelling waves) or evanescent (tunneling) according to whether the term in square brackets is 
positive or negative respectively. A particular plane wave solution of Eq.©, corresponding to an 
incoming wave in region I, is 

z<0 : exp( i n, cos 9 kz) + R exp(— i n 1 cos 8 kz) , 
< z < L : F exp ( i n^h 2 — sin 2 9 kz) + G exp (— i n^h 2 — sin 2 9 kz) , 
z > L : T exp( i n 1 cos 9 kz) , 

with h = rinlrij and i?, F, G and T determined by the boundary conditions. The boundary conditions 
can be set as in quantum mechanics by noting that for piece- wise discontinuities in the "potential" n(z) 
the function U(z) and its first derivative U'(z) must be continuous at the boundaries, or, equivalcntly, 
since the magnetic field is proportional to the derivative of the electric field, by the continuity of these 
fields across the boundaries. Note that the exponentials in region II are oscillatory when n 2 > sin 2 9 
and evanescent for h 2 < sin 2 9. Thus, if n > 1 the solutions will always yield a propagation wave in 
region II of the y-z plane with direction given by Snell's law. When h < 1 both types of solutions 
exist, with tunneling occurring when sin# > n, i.e. for incident angles greater then a critical value 9 C 
(sin 9 C = h). 

The comparison with non-relativistic quantum mechanics also implies that the y-dependence re- 
places the time dependence of the latter. This accounts for the constant term in front of U (z) in 
Eq.©. Another identity is the condition (which we prove in the next section) that 

\R\ +\T\ =1. 

In quantum mechanics this implies conservation of probability [9] while here it implies conservation 
of energy [10]. 
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II. REFLECTION AND TRANSMISSION COEFFICIENTS 

Consider first diffusion (n > sin#) and treat the continuity conditions for U(z) and U'(z) at the two 
interfaces (z = 0, L) independently. Let r and t be the coefficients at the z = interface, then 

r = ^cos 9 — \J h 2 — sin 2 9 J / (cos 9 + \J h 2 — sin 2 9 ^ , 

t = 2 cos 9/ (cos9 + \/n 2 - sin 2 9 ) . (6) 

For a wave travelling from region 11 to region I (e.g. a wave reflected from the z = L interface) the 
corresponding coefficients f Q and t are 

To = — To j 

t = y'n 2 - sin 2 9 t / cos9 . (7) 

Note that for diffusion all the coefficients in Eq.© and ([7|) are real. At the z = L interface, we need 
only consider waves impinging from the left since, for our choice of particular solution, there is no 
incoming wave from the right in region III. Thus, the only reflection and transmission coefficients, are 



r£ = f„ exp ( 2 i n 1 \/ h 2 ~ sin 2 9 kLj 
t L = t exp 



l 71, 



! ( v " 2 — sin 2 ® — cos & ) kL 



(8) 



Now, we may calculate the R and T coefficients by summing individual multiple reflection contribu- 
tions, e.g. the first contribution to R will be r , the second will be i r L anc ^ so f° r th, 

R = r + t rit + < r L f r L i + ... + t r L (f r^) n i + ... (9) 

The series converges because < 9 < ir/2 implies that IrorrJ < I. Summing, we find 

R=r +t r L io/ (l-f r h ) . (10) 

In the same way, 

T = t t L + t f r h t L + ... + t (f r h ) n t L + ... = t t L / (1 - forjj ■ (11) 
Consequently, using the identities r = —r a and r 2 + t t = 1, we find 

\Rf = \r () +r L \ 2 /\l + r a r L \ 2 and \T\" = \t t h \ 2 / |1 + r r h \ 2 . (12) 
It follows after a little algebra that 

\R\ 2 +\Tf = l, (13) 

as anticipated. With the above expressions for R and T, the so called wave limit or total coherence, 
identical to the quantum mechanics results, we reproduce the standard phenomena of resonance when 
\T\ =1. This occurs when the phase in is such that 



= — r n x \J h 2 — sin 2 9 kL = n ir . 

There is however another way to interpret the series expansions for R and T. Let us introduce it by 
simply observing a numerical fact. If we (modulus) square the individual terms in the scries and then 

2 2 

add, we find a different \R\ and |T| , 

\Rn\ 2 = rl + (t r t ) 2 / (1 - <) * \K\ and £ |T n | 2 = (t t a f / (1 - <) ± \T\ 2 , (14) 

n n 

with conservation of energy in the particle limit, where the interference between individual n ampli- 
tudes is null, 

^2(\Rnf +T n \ 2 ) =1. (15) 
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These two limits are easily explained. The former wave limit occurs when all the contribution overlap 
as is the case of plane waves. The second particle limit occurs when no overlapping occurs (see the 
next section) . In quantum mechanics this latter limit corresponds to wave-packets small compared to 
the barrier width (L). The details of how the wave packets are created is not important. It is the limit 
when the time taken for a wave packet to travel back and forth in region II is sufficient to separate 
the individual reflected and/or transmitted wave packets. There arc of course intermediate cases of 
partial overlap. Notice that in the particle limit there is no resonance phenomena. Similarly, for any 

2 

given localized optical transverse electric beam, we can calculate \T\ for various L values (sec next 
section) and see the transition from typical oscillatory (resonance) shape to a constant (particle) limit. 



III. LOCALIZED BEAMS (NUMERICAL ANALYSIS) 

In quantum mechanics the particle limit is obtained by considering narrow (compared to the barrier 
width) wave packets [1,2,11,12]. This is done by integrating the plane wave results with, say, a gaussian 
function in particle momentum. The optical equivalent is to integrate over the incoming angle 9 and 
again this can be done with a gaussian in 9 (see below) . In situations in which tunneling may occur one 
should formally limit the allowed values of 9 to either the diffusion or tunneling regions. In practice a 
strongly peaked dependence around a mean 9 value say 9 = 9 is sufficient as long as Q is sufficiently 
removed from the critical angle 9 C . 

We shall use for our numerical calculations the following gaussian function 



9(a) = ~ — exp 

(27T) 



(a - a ) 5 



(16) 



with a = cos 6? (a = cos9 ) and 5 = n,kd. 

The integration over angles about 9 produces a spatial localization in z and y (the analog of a 
quantum mechanics wave packet). We recall that for our optical study all results are time independent 
(stationary). The localized distributions in z and y are for the incoming, reflected and transmitted 
beams given by 

i 

E^inciy, z) = I da g(a) cxp(i a 5 z d ) cxp(i y/l - a 2 S y d ) , 

J 

1 

E 1 , re f(y,z)= / da R(a) g(a) exp(- iaS z d ) exp(z y/l - a 2 6y d ) , 



E 1} tra{y,z)= / da T(a) g(a) exp(i a S z d ) exp(i a/1 - a 2 8y d ) , 



with y d = y/d, z d = z/d and L d = L/d. 

For diffusion phenomena, and when the beam localizations are smaller than the dimension L of 
region II, wc obtain a zebra- like structure sketched in Fig. 1(b). The various reflected (transmitted) 
beams are separated in the y-z plane. No interference occurs between them and consequently no 
resonance phenomena exists. As in quantum mechanics, these multiple structures do not occur in 
tunneling phenomena [3,4,13], see Fig. 1(a). Indeed for tunneling the sum J2 n \Rn\ diverges, as does 
\Tn\ j thus the individual terms cannot be identified with physical probabilities. In tunneling only 
one reflected and transmitted wave exists. The calculation of R n and T n is then, at best, a technique 
for the determination of R and T. 

We can exhibit these differences graphically by considering a narrow beam incident at an angle 
9 a = 7r/4. Two cases for n will be considered h = 1/2 and h = y/3/2. The choice of a narrow 
beam (5 = 50) guarantees that the gaussian angle distribution, centered in it/ 4, is practically zero 
for 9 < 7r/6 (the critical angle for n = 1/2) and for 9 > 7r/3 (the critical angle for n = y/3/2). 
Consequently, for h = 1/2 we have "tunneling" and for h = V3/2 we have multiple diffusion. In Fig. 2 
and 3, we display, for the diffusion case, plots of {E^ against z for fixed y and against y for fixed z. 
We readily see the multiple beams. Comparison of Fig. 2(a) and 2(b) shows that the reflected beams 
remain separated in z. 
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IV. CONCLUSIONS 

In this paper, we have studied the behavior of a localized beam in a stratified medium. The localiza- 
tion is achieved by integrating over the incidence angle. Depending on the value of n, we have two 
phenomena. One is the formation of multiple beams, the other occurring for tunneling yields a single 
reflected and transmitted beam. We have shown some examples of these phenomena. In both cases 
resonance effects are absent. If we substitute the y-axis with the time axis, we replicate the results of 
multiple diffusion and/or tunneling in non-relativistic quantum mechanics. There are of course some 
significant difference. Foremost, the absence of H and the interpretation of \R\ and |T| in terms of 
energy probabilities. The wave packets in quantum mechanics move in time. In our optical model 
all results are stationary. This is a significantly useful feature, since time measurements are all but 
impractical in quantum mechanics [4] . 

However, the analogy allows us to anticipate some further consequences for localized optical beams. 
For example, while the resonance phenomenon in tunneling is absent for a single barrier, it surprisingly 
reappears [2] for twin or even multiple identical barriers (always in the tunneling regime). Again this 
is a consequence of interference. Thus, if the size of the barriers, including the inter-barrier distances, 
is much larger than the incoming beam, resonance will not occur and multiple beams, caused by 
reflections between barriers, will appear. The ephemeral nature of probability densities makes an 
optical analogy simpler to create and study. It is even conceivable that questions related to tunneling 
times for which there are diverse definitions [4] and the related Hartman effect [3] could be studied 
experimentally. Another potential source of study is the effect of localization on theoretical predictions 
almost always based upon a plane wave analysis. 

The results of this paper, with the wave and particle limits, clearly demonstrate that these effects 
can be significant. In particle physics the effect of wave packets upon phenomena such as neutrino 
oscillations and oscillation phenomena in general, have little or no possibility of experimental testing. 
Perhaps through optics, wc may experiment with some of these questions. 

We have exhibited in our numerical analysis, and shown in our graphs, the absence of a particle 
limit in the case of tunneling. The theoretical method of calculation we have used, based upon the 
sum of individual contributions still works, but only if interpreted as an analytic continuation of the 
diffusion case. The infinite scries in tunneling formally diverges. This is most simply seen by the fact 
in this case r is complex with unitary modulus. 
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Fig. 1. Tunneling (a) and diffusion (b) of a localized optical beam by a dielectric film. The localization is 
achieved by integrating over the incidence angle (narrow gaussian distribution around the angle 8 a — 7r/4). 
For diffusion phenomena, when the beam localizations are smaller than or of the order of the dimension of 



Stefano De Leo, Pietro P. Rotelli: localized beams and dielectric barriers 



7 



- y = 6 d 



0.015 - 



0.010 - 



0.005 - 







- y = 5d 



0.015 
0.010 
0.005 h 




0.015 - 



0.010 - 



0.005 - 







"i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 



2 

\E 1 (y,z)\ vs. z/d 





- y = — 5 d 



J I L 




J I L 



Fig. 2 



(a) - 





(b) - 



(c) " 



nn/ni = V3/2 
O = tt/4 
L = d 
riikd = 50 

J I I I I I I L 



-10 -8 -6 -4 -2 2 4 6 8 10 

2 

Fig. 2. For the diffusion case, the plots of against z for fixed y show the multiple beams and the 

localization in the z-axis. 
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Fig. 3. For the diffusion case, the plots of \E\\ against y for fixed z show the multiple beams and the 
localization in the y-axis. 
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(N : I. INTRODUCTION 

, There exist many unanswered questions in potential theory quantum mechanics. Amongst these is 

' the existence of multiple diffusion phenomena [1,2], the Hartman effect with its apparent violation 

\ of causality [3,4], the importance, if any, of wave packets in oscillation phenomena [5,6]. Most of 

these lack direct experimental measurements. It is therefore extremely instructive to study that class 
of Maxwell equations which are analogous to the Schrodinger equation. These analogies are well 
known in optics e.g. we often find references to tunneling and resonance phenomena. However, their 
5^ ' relevance to quantum mechanics has not been fully exploited. This paper studies an example of the 

above analogy. 

From the Maxwell equations, we can obtain differential equations which the electric and the 
magnetic vector must separately satisfy [7]. For example, for the electric field E, in the case of no 
charges or currents, one has 

v'-E- ^ d tt E+ (Vlnu) x (V x E) + W(E ■ Vine) = . (1) 

c 2 

The corresponding equation for the magnetic field H is obtained by making the changes e fi 
and E — * H . We shall confine our attention to the study of a medium characterized by a real (no 
attenuation) refractive index whose properties are constant throughout each plane perpendicular to 
the chosen direction z, stratified medium [8], 

n(z) = { n, for z < , n u for < z < L , n UI = n 1 for z > L } , (2) 

and for which [i assumes the same value in all three regions. By taking the plane of incidence to be 
the y-z plane, for a monochromatic, exp(— iuot), transverse electric wave (E 2 , 3 = 0), Eq.(fT]) reduces to 

d yy E 1 {y,z) + d zz E 1 (y,z)+n 2 {z)k 2 E 1 {y 1 z) = , (3) 

with k = cu/c and where we have used the fact that V • {eE) = implies that E 1 is a function of y and 
z only. The components of the magnetic vector can be determined by using V x E — — dt (fiH)/c, 

{H 2 (y,z),H 3 (y,z)}^^-{-d z E 1 (y,z),d y E 1 (y,z)} . (4) 
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The geometry of our problem is schematically represented in the following picture: 




z 
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Looking for a separable solution of the form 

E 1 (y, z) = U(z) exp(? n, sin 9 ky) , 

with 9 representing the incidence angle, we obtain the following second-order linear differential equa- 
tion for U(z), 

U"(z) + [n 2 (z)-n 2 +n 2 cos 2 6} k 2 U(z) = . (5) 

This equation is formally identical to the one-dimensional Schrodinger equation when the factor 
which multiplies U(z) is replaced by 2m(E — V)/h . As for Schrodinger the solutions are oscillatory 
(travelling waves) or evanescent (tunneling) according to whether the term in square brackets is 
positive or negative respectively. A particular plane wave solution of Eq.©, corresponding to an 
incoming wave in region I, is 

z<0 : exp( i n, cos 9 kz) + R exp(— i n 1 cos 8 kz) , 
< z < L : F exp ( i n^h 2 — sin 2 9 kz) + G exp (— i n^h 2 — sin 2 9 kz) , 
z > L : T exp( i n 1 cos 9 kz) , 

with h = rinlrij and i?, F, G and T determined by the boundary conditions. The boundary conditions 
can be set as in quantum mechanics by noting that for piece- wise discontinuities in the "potential" n(z) 
the function U(z) and its first derivative U'(z) must be continuous at the boundaries, or, equivalcntly, 
since the magnetic field is proportional to the derivative of the electric field, by the continuity of these 
fields across the boundaries. Note that the exponentials in region II are oscillatory when n 2 > sin 2 9 
and evanescent for h 2 < sin 2 9. Thus, if n > 1 the solutions will always yield a propagation wave in 
region II of the y-z plane with direction given by Snell's law. When h < 1 both types of solutions 
exist, with tunneling occurring when sin# > n, i.e. for incident angles greater then a critical value 9 C 
(sin 9 C = h). 

The comparison with non-relativistic quantum mechanics also implies that the y-dependence re- 
places the time dependence of the latter. This accounts for the constant term in front of U (z) in 
Eq.©. Another identity is the condition (which we prove in the next section) that 

\R\ +\T\ =1. 

In quantum mechanics this implies conservation of probability [9] while here it implies conservation 
of energy [10]. 



Stefano De Leo, Pietro P. Rotelli: localized beams and dielectric barriers 



3 



II. REFLECTION AND TRANSMISSION COEFFICIENTS 

Consider first diffusion (n > sin#) and treat the continuity conditions for U(z) and U'(z) at the two 
interfaces (z = 0, L) independently. Let r and t be the coefficients at the z = interface, then 

r = ^cos 9 — \J h 2 — sin 2 9 J / (cos 9 + \J h 2 — sin 2 9 ^ , 

t = 2 cos 9/ (cos9 + \/n 2 - sin 2 9 ) . (6) 

For a wave travelling from region 11 to region I (e.g. a wave reflected from the z = L interface) the 
corresponding coefficients f Q and t are 

To = — To j 

t = y'n 2 - sin 2 9 t / cos9 . (7) 

Note that for diffusion all the coefficients in Eq.© and ([7|) are real. At the z = L interface, we need 
only consider waves impinging from the left since, for our choice of particular solution, there is no 
incoming wave from the right in region III. Thus, the only reflection and transmission coefficients, are 



r£ = f„ exp ( 2 i n 1 \/ h 2 ~ sin 2 9 kLj 
t L = t exp 



l 71, 



! ( v " 2 — sin 2 ® — cos & ) kL 



(8) 



Now, we may calculate the R and T coefficients by summing individual multiple reflection contribu- 
tions, e.g. the first contribution to R will be r , the second will be i r L anc ^ so f° r th, 

R = r + t rit + < r L f r L i + ... + t r L (f r^) n i + ... (9) 

The series converges because < 9 < ir/2 implies that IrorrJ < I. Summing, we find 

R=r +t r L io/ (l-f r h ) . (10) 

In the same way, 

T = t t L + t f r h t L + ... + t (f r h ) n t L + ... = t t L / (1 - forjj ■ (11) 
Consequently, using the identities r = —r a and r 2 + t t = 1, we find 

\Rf = \r () +r L \ 2 /\l + r a r L \ 2 and \T\" = \t t h \ 2 / |1 + r r h \ 2 . (12) 
It follows after a little algebra that 

\R\ 2 +\Tf = l, (13) 

as anticipated. With the above expressions for R and T, the so called wave limit or total coherence, 
identical to the quantum mechanics results, we reproduce the standard phenomena of resonance when 
\T\ =1. This occurs when the phase in is such that 



= — r n x \J h 2 — sin 2 9 kL = n ir . 

There is however another way to interpret the series expansions for R and T. Let us introduce it by 
simply observing a numerical fact. If we (modulus) square the individual terms in the scries and then 

2 2 

add, we find a different \R\ and |T| , 

\Rn\ 2 = rl + (t r t ) 2 / (1 - <) * \K\ and £ |T n | 2 = (t t a f / (1 - <) ± \T\ 2 , (14) 

n n 

with conservation of energy in the particle limit, where the interference between individual n ampli- 
tudes is null, 

^2(\Rnf +T n \ 2 ) =1. (15) 
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These two limits are easily explained. The former wave limit occurs when all the contribution overlap 
as is the case of plane waves. The second particle limit occurs when no overlapping occurs (see the 
next section) . In quantum mechanics this latter limit corresponds to wave-packets small compared to 
the barrier width (L). The details of how the wave packets are created is not important. It is the limit 
when the time taken for a wave packet to travel back and forth in region II is sufficient to separate 
the individual reflected and/or transmitted wave packets. There arc of course intermediate cases of 
partial overlap. Notice that in the particle limit there is no resonance phenomena. Similarly, for any 

2 

given localized optical transverse electric beam, we can calculate \T\ for various L values (sec next 
section) and see the transition from typical oscillatory (resonance) shape to a constant (particle) limit. 



III. LOCALIZED BEAMS (NUMERICAL ANALYSIS) 

In quantum mechanics the particle limit is obtained by considering narrow (compared to the barrier 
width) wave packets [1,2,11,12]. This is done by integrating the plane wave results with, say, a gaussian 
function in particle momentum. The optical equivalent is to integrate over the incoming angle 9 and 
again this can be done with a gaussian in 9 (see below) . In situations in which tunneling may occur one 
should formally limit the allowed values of 9 to either the diffusion or tunneling regions. In practice a 
strongly peaked dependence around a mean 9 value say 9 = 9 is sufficient as long as Q is sufficiently 
removed from the critical angle 9 C . 

We shall use for our numerical calculations the following gaussian function 



9(a) = ~ — exp 

(27T) 



(a - a ) 5 



(16) 



with a = cos 6? (a = cos9 ) and 5 = n,kd. 

The integration over angles about 9 produces a spatial localization in z and y (the analog of a 
quantum mechanics wave packet). We recall that for our optical study all results are time independent 
(stationary). The localized distributions in z and y are for the incoming, reflected and transmitted 
beams given by 

i 

E^inciy, z) = I da g(a) cxp(i a 5 z d ) cxp(i y/l - a 2 S y d ) , 

J 

1 

E 1 , re f(y,z)= / da R(a) g(a) exp(- iaS z d ) exp(z y/l - a 2 6y d ) , 



E 1} tra{y,z)= / da T(a) g(a) exp(i a S z d ) exp(i a/1 - a 2 8y d ) , 



with y d = y/d, z d = z/d and L d = L/d. 

For diffusion phenomena, and when the beam localizations are smaller than the dimension L of 
region II, wc obtain a zebra- like structure sketched in Fig. 1(b). The various reflected (transmitted) 
beams are separated in the y-z plane. No interference occurs between them and consequently no 
resonance phenomena exists. As in quantum mechanics, these multiple structures do not occur in 
tunneling phenomena [3,4,13], see Fig. 1(a). Indeed for tunneling the sum J2 n \Rn\ diverges, as does 
\Tn\ j thus the individual terms cannot be identified with physical probabilities. In tunneling only 
one reflected and transmitted wave exists. The calculation of R n and T n is then, at best, a technique 
for the determination of R and T. 

We can exhibit these differences graphically by considering a narrow beam incident at an angle 
9 a = 7r/4. Two cases for n will be considered h = 1/2 and h = y/3/2. The choice of a narrow 
beam (5 = 50) guarantees that the gaussian angle distribution, centered in it/ 4, is practically zero 
for 9 < 7r/6 (the critical angle for n = 1/2) and for 9 > 7r/3 (the critical angle for n = y/3/2). 
Consequently, for h = 1/2 we have "tunneling" and for h = V3/2 we have multiple diffusion. In Fig. 2 
and 3, we display, for the diffusion case, plots of {E^ against z for fixed y and against y for fixed z. 
We readily see the multiple beams. Comparison of Fig. 2(a) and 2(b) shows that the reflected beams 
remain separated in z. 
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IV. CONCLUSIONS 

In this paper, we have studied the behavior of a localized beam in a stratified medium. The localiza- 
tion is achieved by integrating over the incidence angle. Depending on the value of n, we have two 
phenomena. One is the formation of multiple beams, the other occurring for tunneling yields a single 
reflected and transmitted beam. We have shown some examples of these phenomena. In both cases 
resonance effects are absent. If we substitute the y-axis with the time axis, we replicate the results of 
multiple diffusion and/or tunneling in non-relativistic quantum mechanics. There are of course some 
significant difference. Foremost, the absence of H and the interpretation of \R\ and |T| in terms of 
energy probabilities. The wave packets in quantum mechanics move in time. In our optical model 
all results are stationary. This is a significantly useful feature, since time measurements are all but 
impractical in quantum mechanics [4] . 

However, the analogy allows us to anticipate some further consequences for localized optical beams. 
For example, while the resonance phenomenon in tunneling is absent for a single barrier, it surprisingly 
reappears [2] for twin or even multiple identical barriers (always in the tunneling regime). Again this 
is a consequence of interference. Thus, if the size of the barriers, including the inter-barrier distances, 
is much larger than the incoming beam, resonance will not occur and multiple beams, caused by 
reflections between barriers, will appear. The ephemeral nature of probability densities makes an 
optical analogy simpler to create and study. It is even conceivable that questions related to tunneling 
times for which there are diverse definitions [4] and the related Hartman effect [3] could be studied 
experimentally. Another potential source of study is the effect of localization on theoretical predictions 
almost always based upon a plane wave analysis. 

The results of this paper, with the wave and particle limits, clearly demonstrate that these effects 
can be significant. In particle physics the effect of wave packets upon phenomena such as neutrino 
oscillations and oscillation phenomena in general, have little or no possibility of experimental testing. 
Perhaps through optics, wc may experiment with some of these questions. 

We have exhibited in our numerical analysis, and shown in our graphs, the absence of a particle 
limit in the case of tunneling. The theoretical method of calculation we have used, based upon the 
sum of individual contributions still works, but only if interpreted as an analytic continuation of the 
diffusion case. The infinite scries in tunneling formally diverges. This is most simply seen by the fact 
in this case r is complex with unitary modulus. 
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Fig. 1. Tunneling (a) and diffusion (b) of a localized optical beam by a dielectric film. The localization is 
achieved by integrating over the incidence angle (narrow gaussian distribution around the angle 8 a — 7r/4). 
For diffusion phenomena, when the beam localizations are smaller than or of the order of the dimension of 
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Fig. 2. For the diffusion case, the plots of against z for fixed y show the multiple beams and the 

localization in the z-axis. 
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Fig. 3. For the diffusion case, the plots of \E\\ against y for fixed z show the multiple beams and the 
localization in the y-axis. 



